Competing orders in Naa;Co02 from strong correlations on a two-particle level 
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Based on dynamical mean-field theory with a continuous-time quantum Monte-Carlo impurity 
solver, static as well as dynamic spin and charge susceptibilites for the phase diagram of the sodium 
cobaltate system Naa;Co02 are discussed. The approach includes important vertex contributions to 
the q-dependent two-particle response functions by means of a local approximation to the irreducible 
vertex function in the particle-hole channel. A single-band Hubbard model suffices to reveal several 
charge- and spin-instability tendencies in accordance with experiment, including the stabilization of 
an effective kagome sublattice close to a;=0.67, without invoking the doping-dependent Na-potential 
landscape. The in-plane antiferromagnetic-to-ferromagnetic crossover is additionally verified by 
means of the computed Korringa ratio. Moreover an intricate high-energy mode in the transverse 
spin susceptiblity is revealed, pointing towards a strong energy dependence of the effective intersite 
exchange. 

PACS numbers: 71.27.+a, 71.30.+h, 71.10.Fd, 75.30.Cr 



The investigation of finite-temperature phase dia- 
grams of realistic strongly correlated systems is a quite 
formidable task due to the often tight competition be- 
tween various low-energy ordering instabilities. In this 
respect the quasi-twodimensional (2D) sodium cobal- 
tate system Na2.Co02 serves as a notably challenging 
case [Il[2]. Here a:S[0,l] nominally mediates between the 
Co4+(3d5,S'=l/2) and Co^+ {icf , S=Q) low-spin states. 
Thus the Na ions provide the electron doping for the 
nearly filled t2g states of the triangular C0O2 layers up 
to the band-insulating limit x—1. Coulomb correlations 
with a Hubbard C/ up to 5 eV for the manifold of 
bandwidth W^~1.5 eV [3] are revealed from photoemis- 
sion Hence with U /W^l the frustrated metallic sys- 
tem is definitely placed in the strongly correlated regime. 

Various different electronic phases and regions for tem- 
perature T vs. doping x are displayed in the exper- 
imental sodium cobaltate phase diagram (see Fig. [ij. 
For instance a superconducting dome (Tc'--^4.5K) sta- 
bilized by intercalation with water close to x=0.3 [5]. 
Pauli-like magnetic susceptibility is found in the range 
x<0.5 T with evidence for 2D antiferromagnetic (AFM) 
correlations i2ii6j. For a;>0.5 spin fluctuations and in- 
creased magnetic response show up for 0.6<a;<0.67, in- 
cluding the evolution to Curie- Weiss (CW) behavior [T] 
for 0.6<a;<0.75, and the eventual onset of in-plane fer- 
romagnetic (FM) order. The ordered magnetic structure 
in the doping range 0.75<x<0.9 with Tn-19-27K [MO] 
is of A-type AFM for the FM C0O2 layers. As the lo- 
cal spin-density approximation (LSDA) is not sufficient 
to account for the AFM-to-FM crossover with x [3], ex- 
plicit many-body approaches are needed |11H13| . 

Several theoretical works have dealt with the influ- 
ence of the sodium arrangements on the electronic prop- 
erties of Naa;Co02, both from the viewpoint of disor- 
dered sodium ions [T2] as well as from orderings for cer- 
tain dopings. jT4Hl6] However, wether such sodium pat- 
terns are due to sole (effective) single-particle potentials 



or mainly originating from many-body effects within the 
C0O2 planes is still a matter of debate [T71 [TB . 

In this letter, we report the fact that a large part of 
the electronic (spin and charge) phase diagram of sodium 
cobaltate may be well described within a Hubbard model 
using realistic dispersions, and without invoking the de- 
tails of the sodium arrangement. Thereby most of the 
observed crossovers and instabilities are truly driven by 
strong correlation effects and cannot be described within 
weak-coupling scenarios. The theoretical study is eluci- 
dating the two-particle correlations in the particle-hole 
channel computed within dynamical mean-field theory 
(DMFT) including vertex contributions (for a review see 
e.g. H^). So far the latter have been neglected in 
cobaltate susceptibilities based on LSDA [31] and 
the fluctuation- exchange approximation |22l 123] . Our 
dynamical lattice susceptibilities allow to reveal details 
of the AFM-to-FM crossover with T and of the intrigu- 
ing charge-ordering tendencies, both in line with recent 
experimental data [2, 24 . Moreover, insight in the (a;, q)- 
dependent spin excitations at flnite frequency is provided. 
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FIG. I. (color online) Left: Sketched Naa;Co02 phase dia- 
gram, based on Ref. i2j Right: M- (top) and K- (bottom) 
point ordering on the triangular lattice. 
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Since we are mainly interested in the a:>0.5 part of 
the phase diagram, the low-energy band dispersion of 
sodium cobaltate is described within an aig-like single- 
band approach, justified from photoemission |25| and 
Compton scattering 26J experiments. We primarily fo- 
cus on the in-plane processes on the effective triangu- 
lar Co lattice with tight-binding parameters up to 3rd 
nearest-neighbor (NN) hopping, i.e., {t,t' ,t")={-202, 35, 
29)meV [57] for the 2D dispersion. Albeit intersite 
Coulomb interactions might play a role |13| , the canonical 
modeling was restricted to an on-site Coulomb interac- 
tion [7=5 eV. Our calculations show that already there- 
from substantial nonlocal correlations originate. The re- 
sulting Hubbard model on the triangular lattice is solved 
within DMFT for the local one-particle Green's func- 
tion G{Ti2)^-{Trc{Ti)c^{T2)) with T„«=T„-T„ and Tr 
being the time-ordering operator. The DMFT prob- 
lem is approached with the continuous-time quantum 
Monte Carlo methodology [28l |29] in its hybridization- 
expansion flavor as implemented in the TRIQS 
package. |30| Additionally we implemented the compu- 
tation of the impurity two-particle Green's function [31J 

G'-^^Ti2,T34,Ti4)=-{Trc\Ti)c{T2)c^{T3)c{T4)) tO addrCSS 

explicit electron- electron correlations. In the approx- 
imation of a purely local particle-hole irreducible ver- 
tex, G*^^^ allows to determine also lattice susceptibili- 
ties, [ini [201 [SD 132] These susceptibilities, e.g. for spin 
(s) and charge (c), written as 

E (^i/i...' (*^' T) + v,i,^,,, (zo., q, T)) , (1) 

where w (y) marks bosonic (fermionic) Matsubara fre- 
quencies, consist of two parts. Namely x^/c denotes 
the conventional (Lindhard-like) term, build up from the 
(renormalized) bubble part, which is mainly capable of 
detecting Fermi-surface driven instabilities close to r=0. 
On the contrary, the second part Vg/c^uu' (the vertex 
term) includes properly the energy dependence of the 
response behaviour due to strong local interactions in 
real space. It proves important for revealing, e.g. mag- 
netic instabilities at finite T due to the resolution of the 
two-particle correlations governed by an implicit inter- 
site exchange J. Note that all numerics take advantage 
of the recently introduced Orthogonal Polynomial repre- 
sentation [31] of one- and two-particle Green's functions 
to provide the needed high accuracy and to eliminate 
artifacts often stemming from truncating the Fourier- 
transformed G*^^^ in Matsubara space. 

Within the first Brillouin zone (BZ) of the triangu- 
lar coordination with lattice constant a the coherent F- 
point instability signals FM order in the case of Xs ^"^^ 
phase separation for Xc- Additionally important are here 
the instabilities at the the K- and M-point. The associ- 
ated orderings give rise to distinct sublattice structures 
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FIG. 2. (color online) Static in-plane charge (top) and spin 
(bottom) susceptibility x(0, q, T) with doping at T=386K. 



in real space (cf. Fig. [I]). The A/-point ordering leads 
to a triangular and a kagome sublattice with lattice con- 
stant acfr=2a, while the /C-point ordering establishes a 
triangular and a honeycomb sublattice with aeff='\/3a, 
respectively. 

We will first discuss the static (Xj/^(w=0, q, T)) re- 
sponse (read off from the zeroth bosonic Matsubara fre- 
quency), directly reflecting the system's susceptibility to 
an order of the (q-resolved) type. The cobaltate intra- 
layer charge-susceptibility Xc(Oiqj^) shows pronounced 
features in q space with doping x (see Fig [2|. Close 
to a;=0.3 our single-band modeling leads to increased 
intensity inside the BZ, pointing towards longer-range 
charge-modulation (e.g. 3x3, etc.) tendencies in real 
space. That Nai/3Co02 is indeed prone to such 120°-like 
instabilities has been experimentally suggested by Qian 
et al. [33]. Towards x=0.5 the susceptibility for short- 
range charge modulation grows in Xcj displaying a diffuse 
high-intensity distribution at the BZ edge with a maxi- 
mum at the K-pomt for a;=0.5. No detailed conclusive 
result on the degree and type of charge ordering for the 
latter composition is known from experiment, however 
chain-like charge disproportionation that breaks the tri- 
angular symmetry is verified [341 135] . The present single- 
site approach cannot stabilize such symmetry-breakings, 
but an pronounced Xc ^-t the iC-point at least inherits 
some stripe-like separation of the two involved sublat- 
tices. Near x=0.67, the Xc maximum has shifted to the 
Af -point, in line with the detection of an effective kagome 
lattice from nuclear magnetic resonance (NMR) experi- 
ments [24]. For even higher doping, this q-dependent 
structuring transmutes into a F-point maximum, point- 
ing towards known phase-separating tendencies. |36j Fig- 
ure [2] also exhibits the x-dependent intra-layer spin sus- 
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FIG. 3. (color online) Temperature dependence of the spin 
susceptibility at the F, M, and A'-point for x=0.55, 0.58, 
0.82. For the latter, values at the yl-point are also included. 
Vertical lines indicate extrapolated transition temperatures 
for F (FM) spin-ordering and A (A-type AFM) spin-ordering 
respectively for 2:=0.82. The inset shows the doping depen- 
dence of the uniform (q=F) in-plane spin susceptibility Xa 
for various T. Note the largely increased magnitude of Xs for 
x>0.75 in this logscale plot. 



ceptibility, starting with strong AFM peaks at a;=0.3 due 
to X-point correlations. With reduced intensity these 
shift to the Af-point at x=0.5, consistent with differ- 
ent types of spin and charge orderings at this doping 
level 35|- For 2;>0.5, Xs(q, T, 0) first develops broad in- 
tensity over the full BZ, before forming a pronounced 
peak at the F-point above Thus the experimen- 

tally observed in-plane AFM-to-FM crossover in the spin 
response is reproduced. 

Lang et al. [2] revealed from the Na NMR that this 
crossover is T-dependent with x, resulting in an energy 
scale T* below which AFM correlations are favored (com- 
pare Fig. [T]). The slope dT* /dx turns out negative, in 
line with the general argument that FM correlations are 
most often favored at elevated T because of the entropy 
gain via increased transverse spin fluctuations. In this 
respect, Fig.[3]shows the (x, T, q) dependence of the com- 
puted Xs- For a;=0.55, 0.58 a maximum in the F-point 
susceptibility is revealed, which has been interpreted by 
Lang et al. [2] as the criterion for a change in the correla- 
tion characteristics, thereby defining the r*-line. While 
the temperature scale exceeds the experimental value in 
the present mean-field formalism, the qualitatively cor- 
rect doping behavior of the r*-line is obtained. 

Beyond the experimental findings our calculations al- 
low to further investigate the nature of the magnetic 
crossover. Fig. |3] reveals that at lower T and x closer 
to x=0.5 the susceptibility at F is ousted by the one at 



M, while Xs ^f ^ is mostly dispensable. The M suscep- 
tibility can be understood due to the proximity of the 
striped order at a;=0.5,[Tl [TTl [37] which is however not 
realized until much lower temperatures. 

The inset of Fig. [3] follows the T-dependent F-point 
susceptibility through a vast doping range. Note the 
subtle resolution around x—Q.b as well as the large exag- 
geration especially for lower temperatures in the exper- 
imentally verified in-plane FM region. The main panel 
of Fig. [3] additionally shows for a;=0.82 the spin suscep- 
tibility at the A-point (i.e., at fc^ = (0,0,l/2) in the BZ), 
which denotes the A-type AFM order. While F and A 
show CW behavior, the extrapolated transition temper- 
ature however is ~7% higher at A than at F, verifying 
the experimental findings of A-type order [7HTU] . In the 
temperature scan we additionally introduced a nearest- 
layer inter-plane hopping ij^=13 meV [SI [T31 [35]; how- 
ever the previous in-plane results are qualitatively not 
affected by this model extension. Due to known charge 
disproportionation the inclusion of long-range Coulomb 
interactions, e.g., via an inter-site V jl3l 116] . seems rea- 
sonable. This was abandoned in the present single-site 
DMFT approach, resulting generally in reduced charge 
response. Without V , charge fluctuations are substan- 
tially suppressed for large U /W , while the inter-site spin 
fluctuations are still strong due to superexchange. Aside 
from the static response, our method allows also access 
to the dynamic regime. Figure |4] shows the dynami- 
cal transverse spin susceptibilitiy for selected x. Note 
the broad q-dependence and small excitation energy in 
the low-doping regime. In contrast, the FM correlations 
near x=0.82 are reflected by strong paramagnon-like gap- 
less excitation at F combined with very little weight and 
rather high excitation energies at AFM wave- vectors. In- 
terestingly, a comparably strong and sharp iiT-type high- 
energy excitation (~1 eV )for larger x below the onset of 
in-plane FM order is revealed. Its amplitude is strongest 
at a;=0.67 while its energy increases with x and its worth- 
while to note that the mode is neither visible when ne- 




FIG. 4. (color online) Imaginary part of the dynamical spin- 
susceptibility for selected dopings for selected dopings. 
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FIG. 5. (color online) Korringa ratio versus doping for r=580 
K. The experimental data, extracted from Ref. HO] and [2] 
was obtained for lower temperatures. The inset shows the 
evolution of the bubble-diagram contribution from the non- 
interacting (!7=0) to the fully-interacting {U—5 eV) case. 



glecting vertex contributions, nor in a plain triangular 
Hubbard model with NN hopping only. Thus it reflects a 
strong energy dependence of the inter-site exchange cou- 
pling J= J(a;, q, oj), that obviously changes character for 
with q and w. The predicted high-energy feature 
could be probed experimentally and also studied in time- 
dependent measurements. We propose the use of modern 
laser-pulse techniques [31] to address this problem. 

Experimentally, the evidence for significant q^O fluc- 
tuations is drawn [2 HD] from the Korringa ratio 
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where l/Ti is the nuclear relaxation rate, Ks the NMR 
field shift, 7e {in) the electronic (nuclear) gyromagnetic 
ratio and A(q) the hyperfine coupling. Roughly speak- 
ing, 1C>1 signals AFM correlations, /C<1 points to FM 
tendencies in Xs ^'^'^ the term "Korringa behavior" gen- 
erally denotes the regime /C(T)~1. In single-atom unit 
cells, ^(q) becomes q- independent. 

Note that especially l/Ti is numerically expensive, as 
it requires to calculate on many Matsubara frequen- 
cies with subsequent analytical continuation to the real 
frequency axis for contributions beyond the bubble di- 
agram. Figure [5] finally shows the AFM-to-FM correla- 
tion crossover captured by the Korringa-ratio over a wide 
doping range. The overall agreement with experiment is 
conclusive. Relevant deviations in the low-doping regime 
probably originate from the smaller temperatures studied 
in experiment. The difference at a;=0.58 might be of the 
same origin, but since charge ordering occours for a::>0.5 
which was not included explicitly here, neglecting the q- 
dependence of A(q) might be also questionable. |45) One 



can see that the bubble-only calculation yields a nearly 
fiat Korringa ratio with doping, thus fails completely in 
explaining the experimental findings. In particular it 
does not refiect the strong FM correlations for high dop- 
ing. This further proves the importance of strong correla- 
tions on the two-particle level, asking for substantial ver- 
tex contributions. [33] Note that the recently suggested 
lower-energy effective kagome model [H] including the 
affect of charge ordering is not contradicting the present 
modeling. Since here the effective kagome lattice natu- 
rally shows up and also the key properties of the spin 
degrees of freedom seem well described on the original 
triangular lattice. 

In summary, the DMFT computation of two-particle 
observables including vertex contributions based on a re- 
alistic single-band Hubbard modeling for Naj;Co02 leads 
to a faithful phase-diagram examination at larger x, 
including the kagome-like charge-ordering tendency for 
a;~0.67 and the in-plane AFM-to-FM crossover associ- 
ated with a temperature scale T* . Thus it appears that 
many generic cobaltate features are already governed by 
a canonical correlated model, without invoking the de- 
tails of the doping-dependent sodium-potential landscape 
or the inclusion of multi-orbital processes. Of course, 
future work has to concentrate on quantifying further 
details of the various competing instabilities (and their 
mutual couplings) within extended model considerations. 
Beyond equilibrium physics, we predict a strong energy 
dependence of the effective inter-site exchange resulting 
in an if-type high-energy mode around x=0.67, which 
could be probed in experimental studies. 
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